Capillary Waves in Binder's Approach to the 

Interface Tension 



(N . U.-J. Wiese" 

a^ 

^! Universitat Bern, Sidlerstrasse 5, 3012 Bern, Switzerland 

^ ■ February 1, 2008 

> 
O 

o 

a^ 

^ I In Binder's approach the reduced interface tension a of the Ising 

0\ I model in the broken phase is determined from the finite volume effects 

of the partition function Z{M) at fixed total magnetization M. For 
small \M\ the partition function of a system of size L'^ with periodic 
Qh! boundary conditions is dominated by configurations with two interfaces, 

r-| ! such that Z{M) oc exp(— 2crL^^^). Capillary wave fluctuations of the 

interfaces correct this result to Z{M) oc exp{—2crL'^^^) with x = — 1. 
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Abstract 



The knowledge of the pre-exponential behavior allows an improved fit of 
numerical data, and a determination of the interface stiffness. 
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Let us consider the d-dimensional Ising model on a periodic L lattice in the low 
temperature broken phase. Recent numerical simulations of the so-called multimag- 
netical ensemble in the d = 2 and in the d = 3 case have demonstrated that 
Binder's approach to the interface tension is numerically feasible. The numerical 
method computes the partition function as a function of the total magnetization 
M = J^x'^ix), where (t{x) = ±1 are the spin variables. The partition function at 
fixed total magnetization is given by 

a^' a(x)=±l 

where Sm,n is the Kronecker 5-function, P is the inverse temperature and Hla] = 
— J2{xy) '^{^)(^{y) is the Hamilton function with nearest neighbor interactions. When 
\M\ is small, typical configurations consist of two domains (of + and — phase) 
separated by two interfaces, which are closed via periodic boundary conditions. The 
dynamics of the interfaces can be described in a capillary wave model. In the d = 2 
case the interfaces are lines, which one can parametrize as Xi{t) and X2{t) with 
t E [0,L]. A typical configuration is depicted in fig.l. 



Figure 1: A typical configuration consisting of two domains of + and — phase 
separated by two interfaces parametrized as Xi{t) and X2{t). 
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The energy of the configuration can be written as 

'H[xuX2] = £ dt[fL + a{ipl)^JlTi^ + a{ip2)^/lT^], (2) 

where / is the bulk free energy density, (pi = atan(xj) (with Xj = dxi/dt) are the 
angles between the interfaces and the lattice t-axis, and a{(f ) is the angle dependent 
interface tension, which is known analytically . The minimal energy configuration 
contains two fiat parallel interfaces (with (pi = 0) and it has Ti.[xi, X2] = /L^+2a(0)L 
independent of M. The leading finite volume effect of Z{M) is therefore given by 
Binder's result 

Z{M) oc exp(-/3/L2) exp(-2aL), (3) 

where a = Pa{0) is the reduced interface tension. Eq.(0) receives corrections due 
to interface fluctuations. Binder has parametrized this effect by a pre-exponential 
factor L^. One goal of this paper is to compute the value of x. Since we are only 
interested in the next to leading flnite volume effect of Z{M), we may assume that 
the interfaces are almost flat with only small capillary wave fluctuations [^, such 
that ^ 1. Then 

/3n[xi, X2] = dt[f3fL + 2a + ^x? + ^xj], (4) 

where k = P{a{0) + «"(0)) is the reduced interface stiffness. (Note that a'{0) = 0.) 
The area A = Jq dt{xi — X2) between the two interfaces determines the value of 
the total magnetization M = m{2A — L?) of the conflguration, where m is the 
magnetization density. For example, when the enclosed area is maximal {A = L^) 
the magnetization is M = mL^, while for A = one has M = —mL?. For the 
partition function at flxed magnetization we now write the following path integral 

Z{M) = 2 / Vxi{t)Vx2it)6{ 2 / dt{xi - X2) + L^) exp(-/?7^[xl, X2]). (5) 

J m Jo 

The integration is over all pathes (interface locations) xi{t) and X2{t), which are peri- 
odic in t and which do not cross each other. The non-crossing constraint is a natural 
topological consequence, due to the fact that the Xi{t) are interfaces separating the 
+ from the — bulk phase. The factor 2 in front arises because the interfaces can wrap 
around the lattice in two different directions. The system of eq.(§) corresponds to 
two nonrelativistic particles (with coordinates Xi{t) and X2{t)) on a circle of radius 
L, interacting repulsively at short distances. In addition the area enclosed between 
the pathes of the two particles is constraint to A = (M/m+L^)/2. The path integral 
eq.(^) is equivalent to a two-particle Schrodinger equation. To solve the Schrodinger 
equation it is useful to introduce a center of mass coordinate X = (xi + X2)/2 and 
a relative coordinate x = xi — X2. The center of mass motion does not feel the 
constraint. It is described by plane waves with total momenta P = 2ttI/L {I G Z) 
and with energies Ei = (3fL + 2a + P'^/Ak = /SfL + 2a + tiH'^ / kL'^ . The partition 
function factorizes 

Z{M) = 2 exp{-/3fL'^) exp(-2aL) ^ exp{-n^l^ / kL)z{M). (6) 

/GZ 
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Since we are interested in the large volume limit we may replace the sum by an inte- 



gral Z^iGZ 6xp(— vt^/^/kL) ~ (i/exp(— vr^/^/ftL) = y^/tL/vr. The relative motion 
is described by 

z{M) = I Vx{t)5{— - 2 dtx + L^) exp(- T dt^x^) 
J m Jo Jo A 

= _/ deexpHe—) Vx(t)exp(- dt[-x^ + 2ie(x - -)]). (7) 
27r J-oo m J Jo A 2 

In the last step the constraint of fixed magnetization has been rewritten as an 
integral over an additional parameter 6. The path integral for the relative motion 
corresponds to a Schrodinger equation for a wave function \E'„(a;) with an imaginary 
linear potential 

- -^^^^ + W{x - ^)^n{x) = e^{9)^r^{x). (8) 

The energies e„(^) are in general complex. The non-crossing constraint of the inter- 
faces refiects itself in the boundary conditions \E'„(0) = \l/n(-^) = 0. (One may also 
impose other boundary conditions without changing the final result.) The partition 
function of the relative motion is given by 

z{M) = — / deexv{ie—)Y,exv{-tn{9)L), (9) 

27r j-oo m „ 

and the problem reduces to the computation of the ^-dependent energy spectrum. 
Since we are interested in small |M|, we need e„(6') only for large values of \9\. The 
Schrodinger equation (|^) is solved in terms of Airy functions '^n{x) = aAi{zn{x)) + 
bBi{znix)) with z„(x) = {K/Ae^/^itniO) ~ 2ie{x - L/2)). The boundary conditions 
result in the quantization condition Ai(2;„(0))Bi(z„(L)) = Ai(2;„(L))Bi(2;„(0)). For 
large values of |^| the solutions come in pairs with complex conjugate energies and 
one finds 

e„(^) ~ ±i\e\L + {Ae^/Ky/^rnexp{±2m/3), (10) 

where the r„ < are the zeros of the Airy function Ai(r„) = 0. In the large volume 
limit one can use the asymptotic form r„ ~ — (37rn/2)^/^ and one can replace the 
sum over n in eq.(^ by an integral. After some manipulations one finds 

1 r°° M r°° 

z{M) = -^-^ decos{e—) dy^exp{-y/2KL)x 
h^kL'^ Jo m Jo 

[cos{9L'^) cos{V3y/2KL) + sm{eL'^) sm{V3y/2KL)]/9, (11) 

where y = {StttikOL'^)'^/'^. Using dyy/yexp{-y /2kL) cos{\/3y /2kL) = and 
dy^ exY>{—y /2kL) sm{^/?)y /2kL) = {kL)'^/'^ ^ /2 one obtains 

1 [kLI .„ ,M 



z(M) = —J — / decos(e—)sm(eL^)/e 

2L^ \ 71 71 Jo m 




for \M\ < mL'^, and z{M) = for \M\ > mL"^ . Still, in the small \M\ region, 
the partition function is independent of the magnetization. This is because the 
fluctuations of the interfaces are small and do not induce interactions between the 
two interfaces. In total one finds 

Z{M) = exp{-/3fL'^)^exp{-2aL) for |M| < mL^. (13) 

2tcL 

This corresponds to x = —1 in Binder's pre-exponential factor L^. Eq.(]T3|) is reliable 
only for small values of |M|, because only in this region the typical configurations 
consist of two domains separated by two interfaces. When \M\ increases and one 
of the domains shrinks in favor of the other one, the lowest energy configuration 
contains a droplet of one phase in a background of the other phase. The equilibrium 
shape of the droplet follows from the so-called Wulff construction 0. Using the 
analytic expression for the angle dependent interface tension a{ip), this construction 
has been carried out for the d = 2 Ising model in ref. 0. The equilibrium droplet has 
a spherical shape only close to the second order phase transition, where a continuum 
limit is approached and a{(p) becomes angle independent. In this case the region 
where a droplet costs less energy than two parallel interfaces is at \M\ > mL^(l — 
2 /it). In the zero temperature limit, on the other hand, the equilibrium droplet 
becomes quadratic. Then a droplet is energetically favorable for |M| > mL^/2. The 
general case hes between the two extremes. Of course, also a droplet has capillary 
wave fluctuations. They could be handled in the same way as before. However, away 
from the critical region the angle dependence of a{(f ) complicates the situation. In 
any case, we are most interested in the small \M\ region where the results from 
above apply. 

In the numerical simulations of the multimagnetical ensemble one measures 
Z{M) only up to a global normalization. In practice the reduced interface ten- 
sion a is extracted from the ratio of the partition function at the minimum Z{0) 
and at the maximum Z(mL2). At the maximum typical configurations consist of 
just one domain (of either + or — phase) and no interfaces are present. Such pure 
phase configurations have so far not been included in the calculation. It is trivial 
to include these configurations, because they only contribute a (5-function to the 
M-dependent partition function. The 5-function is normalized by the Boltzmann 
factor of the free energy, and we can write 

Z(M) = expi-PfL'M- - L') + 5{- + L') + ^ exp(-2aL)]. (14) 

Of course, also a pure phase has fluctuations in the total magnetization, which have 
not yet been taken into account. It is well known that the (5-functions are smeared 
to Gaussian distributions D{M'/m) = -i/l/27rxL^ exp(— M'^/2m^xL^) in a finite 



volume (see e.g. [^). Here x is the magnetic susceptibility and M' = M ± mL^. 

The same fluctuations should be included in the two phase configurations at small 

\M\. This is achieved by the convolution 

poo M' M> 
ZiM) ^ / d — ZiM-M')Di — ) 
j-oo m m 
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= exp(-/3/L2)[D(--L2) + D(- + L2) + -^exp(-2aL)]. (15) 

This is the final result, which is valid in the two interface regime at small \M\ as 
well as near the pure phase regions M = ^mL"^. At intermediate |M|, in the region 
where droplets dominate, it is not valid and we will not use it there. Most interesting 
for applications to numerical simulations is the ratio 

Z{0) _ fX^exp{-2aL) ford = 2, (16) 



Z{mL^) V 277 

which has an L-independent pre-exponential factor. The L-independence was al- 
ready observed empirically in ref.[|l| by fitting the numerical data. Eq.(|T^) allows 
to determine the interface stiffness k from a fit of numerical data, when x is known. 

Now let us turn to the d = 3 case. It is important to note that 3-dimensional 
interfaces living on a lattice have a roughening transition (see e.g. 0). At very low 
temperatures the interfaces are rigid, with small steplike excitations following the 
lattice structure. At temperatures above the roughening transition, which is well 
below the bulk critical temperature, interfaces become rough and fluctuate more 
freely. Q In particular, the capillary waves arise as soft modes of interface fluctua- 
tions. Here we restrict ourselves to temperatures above the roughening transition. 
First, let us consider a volume L x L x L' with L' -C L. Then capillary waves in the 
short direction cost much energy and are therefore suppressed, such that the prob- 
lem remains quasi 2-dimensional. Repeating the whole calculation with the extra 
factor of L' one obtains 

Z{M) = exp{-(3fL^L')\D{— - L^L') + D{— + L^L') + ^ exp{-2aLL')], (17) 



where now D{M'/m) = ^1/2ttxL'^L' exp{—M''^ /2m^xL'^L'). Again this corre- 
sponds to X = —1 in Binder's pre-exponential factor. When L' increases and the 
problem becomes truely 3-dimensional, eq.(p!7D is no longer correct, because capil- 
lary waves can also arise in the third direction. As a consequence, numerical factors 
will change, but the general form of the L-dependence should remain the same. In 
the L' — i> L limit I therefore conjecture for the partition function ratio 

^^^^ oc VLexp(-2aL2) for d = 3, (18) 



Z(mL 

such that the pre-exponential factor is now L-dependent. This behavior is in good 
agreement with the numerical data of ref. for large volumes. The knowledge 
of the pre-exponential behavior flxes one free parameter of the flt in ref. [||], and 
therefore improves the analysis of the numerical data. In dimensions c? > 4 it is 
generally believed that interfaces are always rigid. Therefore one must be careful 
in generalizing the results to higher dimensions. It is straightforward to generalize 

"'^Note that in two dimensions interfaces are always rough, except at zero temperature. 
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the results to other models. Especially, the form of the L-dependence should be 
universal. 

I would like to thank P. Hasenfratz for very interesting discussions. 1 am also 
indepted to V. Privman for explaining the concept of capillary waves, and for very 
useful remarks. 
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